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The ability to perform complex as well as simple cognitive tasks engages a network 
of brain regions that is mediated by the white matter fiber bundles connecting them. 
Different cognitive tasks employ distinctive white matter fiber bundles. The temporal 
lobe and its projections subserve a variety of key functions known to deteriorate during 
aging. In a cohort of 52 healthy subjects (ages 25-82 years), we performed voxel-wise 
regression analysis correlating performance in higher-order cognitive domains (executive 
function, information processing speed, and memory) with white matter integrity, as 
measured by diffusion tensor imaging (DTI) fiber tracking in the temporal lobe projections 
[uncinate fasciculus (UF), fornix, cingulum, inferior longitudinal fasciculus (ILF), and 
superior longitudinal fasciculus (SLF)]. The fiber tracts were spatially registered and 
statistical parametric maps were produced to spatially localize the significant correlations. 
Results showed that performance in the executive function domain is correlated with DTI 
parameters in the left SLF and right UF; performance in the information processing speed 
domain is correlated with fractional anisotropy (FA) in the left cingulum, left fornix, right 
and left ILF and SLF; and the memory domain shows significant correlations with DTI 
parameters in the right fornix, right cingulum, left ILF, left SLF and right UF. These findings 
suggest that DTI tractography enables anatomical definition of region of interest (ROD 
for correlation of behavioral parameters with diffusion indices, and functionality can be 
correlated with white matter integrity. 
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The ability to perform complex as well as simple cognitive tasks 
engages a network of brain regions mediated by the white matter 
fiber bundles connecting them. Different cognitive tasks employ 
distinctive white matter fiber bundles. For example, memory 
involves the temporal lobe white matter projections connecting 
the hippocampus with other brain structures. 

Aging is accompanied by deterioration of cognitive functions, 
as well as changes in gray and white matter integrity. Cross sec- 
tional and longitudinal behavioral studies have found significant 
declines in various abilities, including encoding and retrieval 
of new memories, working memory, executive functions, pro- 
cessing speed and spatial ability (Craik et al., 1987; Petersen 
et al, 1992; Youngjohn and Crook, 1993; Small, 2001). Several 
imaging studies tried to relate diffusion parameters measured in 
white matter fiber systems and age related cognitive decline. It 
was shown that reduced integrity of specific fiber tracts are dif- 
ferentially correlated with declines in components of executive 
functions, including working memory, problem-solving, categor- 
ical task switching and Stroop color-word interference (Sullivan 
et al., 2006; Madden et al, 2009; Zahr et al., 2009). In cortical 
association fiber tracts, correlation was found with decrement in 
set-shifting performance observed with age (Perry et al., 2009). 



Motor scores were correlated with FA or diffusivity in several lat- 
eral fiber bundles and in several callosal sectors (Sullivan et al., 
2010). DTI parameters in the ILF and IFOF were correlated 
differently with scores in different visual tasks, and posterior pro- 
jections of the corpus callosum (CC) were found to correlate to 
memory and executive functions (Voineskos et al., 2012). Voxel- 
wise correlation yielded a significant negative association between 
FA in the body of the CC and performance in a complicated 
motor task (Koch et al., 2012). 

Most of the above mentioned DTI studies used averaged 
DTI parameters along the fiber tracts. However, DTI param- 
eters can be heterogeneous along the tract and the spatial 
information obtained from the fiber tracking procedure may 
be important for more accurate understanding of the relation- 
ship between cognitive function and white matter integrity. 
Deterioration of brain tissue during aging involves the tempo- 
ral lobe, along with the cognitive functions it subserves. The 
temporal lobe has several white matter projections: the cingu- 
lum is known to be involved in many brain functions, such as 
working memory (Sepulcre et al., 2009), and attention (Nestor 
et al., 2007; Hamilton et al., 2008); The fornix is important for 
learning and formation of new memories (Spiers et al., 2001); 
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The superior longitudinal fasciculus (SLF) has an important 
role in higher brain functions, more particularly language and 
language disorders (Lichtheim, 1885; Damasio and Damasio, 
1980; Tanabe et al, 1987; Catani et al, 2003; Wise, 2003; 
Geldmacher et al., 2007); the inferior longitudinal fasciculus 
(ILF) is known to play an important role in visual memory as 
demonstrated by post mortem studies (Bauer and Trobe, 1984), 
lesion studies (Shinoura et al., 2007) and imaging studies in 
congenital prosopagnosia (Thomas et al., 2009) and face recog- 
nition (Tavor et al., 2010); the uncinate fasciculus (UF) plays 
an important role in the formation and retrieval of memories 
as shown by lesion (Levine et al., 1998) and animal studies 
(Squire et al., 2004). 

In the present paper we used the inter-subject variability 
in various cognitive domains to relate cognitive performance 
and white matter integrity in five temporal projections: the UF, 
fornix, cingulum, ILF, and SLF. We computed voxel-wise cor- 
relations between cognitive performance in higher-order cogni- 
tive domains (executive function, information processing speed, 
and memory) and white matter integrity, as measured by ADC, 
and FA along the different temporal fiber tracts, correcting for 
age and motor ability. 

METHODS 
SUBJECTS 

Subjects were 52 healthy volunteers, age 25-82 years, all right 
handed, 20 males and 32 females. Mean/median age was 51 years. 
Age distribution was approximately uniform across the age range. 
The local Tel Aviv-Sourasky Medical Center Helsinki Committee 
approved the research protocol, and all participants signed an 
informed consent. Subjects had no history of neurological dis- 
eases or psychological disorders, did not use neuropsychiatric 
medication, had no history of drug or alcohol abuse, and had 
intact vision. All subjects were right handed and had at least 
15 years of education. Volunteers with common age-related dis- 
eases such as diabetes, high blood pressure, and high levels of 
cholesterol were not excluded from the cohort. 

NEUROPSYCHOLOGICAL TESTS 

Participants underwent a series of computerized cognitive tests 
(Mindstreams, NeuroTrax Corp., Bellaire, TX) that evaluate per- 
formance across an array of cognitive domains known to dete- 
riorate during aging (memory, executive function, visual spatial 
processing, verbal function, attention, information processing 
speed, and motor skills) and provide measurements of accuracy 
and reaction time (RT) shown to be valid and reliable (Russell, 

2011) . These tests have been used in clinical settings, as well 
as in studies of normal aging relating these neuropsychological 
measurements to, e.g., genetic findings (Thaler et al., 2012), psy- 
chiatric and neurodegenerative disorders (Doniger et al., 2006; 
Plotnik et al., 2011), and imaging parameters (Sasson et al., 
2010, 2012). In the present study, we performed factor analy- 
sis of 44 cognitive outcome parameters. A detailed description 
of the cognitive tests used in the correlation analysis with DTI 
parameters appears below (in brief). A more detailed descrip- 
tion of the cognitive tests was published previously (Sasson et al., 

2012) . 



FACTOR ANALYSIS PROCEDURE 

Factor analysis was performed using SPSS software (SPSS Inc., 
Chicago, IL, USA), as follows: 

1. Factor analysis of 44 cognitive outcome parameter scores 
using principal component analysis (PCA) with orthogonal 
( VARIMAX) rotation and limiting the number of factors to 5 
(although 10 factors obtained eigenvalues higher than 1) that 
explained 65% of the items total variance. The number of 
factors selected corresponded to the point of inflexion in the 
respective scree plot. The percentage of variance explained 
derives from that point. 

2. Based on the rotated component matrix, and as we aimed 
to extract "pure" factors, we dismissed items with a differ- 
ence below 0.1 between the 2 highest loadings. 14 scores were 
dismissed. 

3. Each factor contained cognitive scores with loadings greater 
than 0.4 as accepted in the PCA literature (Stevens, 2002). 

4. Cronbach's alpha was used to assess the internal reliability 
of the factors. Factors with a Cronbach's alpha of <0.8 were 
excluded, resulting in the exclusion of two factors and the 
retention of three. 

5. The three factors retained were comprised of 18 cogni- 
tive scores. The scores and their factor loadings are shown 
in Table 1. Factor loadings that correspond to the specific 
domain are shown in bold. Each factor was assigned a 
name describing its constituent scores: executive function, 
memory, and information processing speed. 



Table 1 | Factor analysis weights of final 20 cognitive measures, 
assigned to each cognitive domain. 





Executive 
function 


Memory 


Information 

processing 

speed 


Stroop Wl RT 


0.82 


-0.15 


0.32 


Stroop Wl RD 


0.80 


-0.14 


-0.22 


GoNoGo SD 


0.80 


-0.17 


0.12 


Stroop NO 11 RT 


0.74 


-0.33 


0.10 


GoNoGo RT 


0.72 


-0.11 


0.33 


Stroop NOI1 SD 


0.62 


0.07 


0.13 


NV Memory im rep3 


-0.20 


0.87 


-0.14 


NV Memory im rep4 


-0.27 


0.86 


-0.18 


NV Memory Delayed 


-0.28 


0.85 


-0.14 


NV Memory im rep2 


-0.21 


0.82 


-0.17 


NV Memory im rep1 


0.037 


0.71 


-0.16 


Verbal Memory im rep 1 


-0.11 


0.47 


-0.33 


Verbal Memory Delayed 


-0.27 


0.45 


-0.27 


lnfoproc2_1 RT 


0.03 


-0.22 


0.84 


Info proc 2_2 RT 


0.15 


-0.16 


0.81 


Info proc 2_3 RT 


0.25 


-0.34 


0.64 


Info proc 2_3 SD 


0.0005 


-0.05 


0.48 



Note: WO/ 7, no interference 7 (letter color); Wl, with interference; RT, reaction 
time; RD, response time difference, Wl minus NOI2 (word meaning); SD, stan- 
dard deviation; Acc, accuracy; NV, non-verbal; rep, repetition; im, immediate; Info 
proc X_Y, information processing X, number of digits, Y, rate of presentation. 
Bold numbers indicate factor loadings corresponding to each specific domain. 
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6. Calculation of factor scores: for each factor, the z-scores of 
the tests comprising it were averaged to give the factor score. 
These factor scores served as a covariate input for the DTI 
indices. 

DESCRIPTION OF THE COGNITIVE TESTS COMPRISING THE COGNITIVE 
FACTORS 

(a) The executive function factor included scores from the Go- 
NoGo response inhibition test and Stroop interference test. 
The Go-NoGo test is a continuous performance test of sim- 
ple RT and response inhibition. In the Stroop interference test 
performed here, the conflicting information provided by the 
meaning of a word and the color of its letters leads to a decre- 
ment in performance relative to the test phases where there 
is no conflict. In total, the test includes three phases — phase I 
(no interference: letter color), phase II (no interference: word 
meaning), and phase III (interference between letter color 
and word meaning). The executive function domain is an 
average of RT in phase III (interference), response time differ- 
ence (RD) between RT in phase III and phase II, RT in phase 
I (no interference: letter color) and standard deviations (SD) 
in phase I (no interference: letter color). 

(b) The memory domain included non-verbal and verbal mem- 
ory tests. The verbal memory test measures immediate and 
delayed recognition memory for verbal paired associates. 
In the non-verbal memory test, participants are presented 
with eight black-and-white drawn pictures of simple com- 
mon objects or shapes and are instructed to remember their 
orientation. Each of the tests included four repetitions of 
the study phase and a recognition test. This is followed by 
a delayed recognition test after two intervening tests. The 
memory domain included accuracy on the four immediate 
repetitions and delayed phase of the non-verbal memory test 
and accuracy on the first immediate repetition and delayed 
phase of the verbal memory test. 

(c) The staged information processing speed test measures infor- 
mation processing at increasing levels of complexity. The test 
is comprised of three blocks of information processing load: 
single digits, two-digit arithmetic problems (e.g., 5 — 1), and 
three-digit arithmetic problems (e.g., 3 + 2 — 1). For each 
of these three blocks, stimuli are presented at three differ- 
ent rates (speed levels), incrementally increasing as testing 
continues. The information processing speed domain is an 
average of the RT in block 2 (two-digit arithmetic problems) 
at the three different speed levels and SD of the RT in block 
2, level 3. 

MRI ACQUISITION 

MR imaging was performed on a 3T (GE) MRI system at the Tel 
Aviv-Sourasky Medical Center. The MRI protocol included con- 
ventional anatomy sequences and DTI acquired with a standard 
head-coil. 

Conventional anatomy sequences 

Ti -weighted images: 3D spoiled gradient recalled echo (SPGR) 
sequence with the following parameters: 66 axial slices, TR/TE = 
400/3.2 ms, and resolution of 1 x 1 x 2 mm 3 , scan time of 4 min. 



DTI protocol 

Spin-echo diffusion weighted echo planar imaging (DW-EPI) 
sequence was performed with 48 axial slices and resolution 
of 2.5 x 2.5 x 2.5 mm 3 . Diffusion parameters were: TR/TE = 
25/19 ms, fo-value of 1000 s/mm 2 acquired with 19 gradient direc- 
tions. The DW-EPI sequence was gated to the cardiac cycle with 
TR of 30 R-R intervals and TE of 88 ms. 

IMAGE ANALYSIS 

Correction of head motion image artifacts, normalization, and 
statistical analysis were performed using the SPM software (ver- 
sion 2, UCL, London, UK). The image analysis routine of the 
diffusion tensor data included the following steps: 

1. All of the DWI images were coregistered using SPM2 (UCL, 
London, UK) to correct for head motion. Gradient orienta- 
tions were compensated prior to the b-matrices calculation 
to account for the rotation component of the registration. 

2. Spatial normalization was applied: images were normalized 
to the MNI coordinates, using non-linear deformations in 
SPM2, and compensating the gradients orientations for 
the rotation component of the affine-transformation 
that is the closest to the non-linear deformations. 
Normalization was performed using a 12-parameter 
affine transformation followed by non-linear transforma- 
tions, with 1176 parameters describing each deformation 
field. 

FIBER TRACKING PROCEDURE 

Tensor fitting was performed using the free water elimination 
method (Pasternak et al., 2009). The free water elimination 
method allows correcting for water contamination, especially in 
brain tissue adjacent to the lateral ventricles, for example the 
fornix. It is achieved by fitting a bi-tensor model for which a 
mathematical framework is introduced to stabilize the fitting. In 
this routine, we eliminate a tensor component that fits the diffu- 
sion properties of free water and use for tracking only the fiber 
component of the tensor. The free water elimination is important 
when focusing on elderly subjects as in the present work, since 
atrophy and enlarged ventricles might result in lower FA values 
in white matter structures such as the fornix, which is adjacent to 
the lateral ventricles. 

The tensors obtained were spectrally decomposed to their 
eigen-components. The eigenvalues were used to calculate FA 
maps (Basser and Pierpaoli, 1998). Tractography was applied 
using the principal eigenvectors and FA: the brute force FACT 
algorithm was used to generate the fiber coordinates, terminat- 
ing at voxels with FA lower than 0.2 or following tract orientation 
change higher than 60°. Fibers that passed through a manually 
chosen seed region of interest (ROI) were plotted. The fibers 
were plotted as streamlines using Matlab (Mathworks, Natick, 
MA). Once a subset of fibers had been found, a visitation map 
was generated, indicating voxels with at least one streamline 
passing through. The masks obtained were overlaid over the Bo 
image. Overall 10 fibers tracts were plotted for each subject: five 
temporal projections (UF, fornix, cingulum, SLF, ILF), in both 
hemispheres. 
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Representative seed ROIs are shown in Figure 1, as well as the 
reconstructed fiber. The seed ROIs were chosen based on visual 
comparison with the MRI atlas of human white matter (Mori and 
Crain, 2005). The seed ROIs were drawn on the DTI color map 
for each temporal projection in the following manner: 

1. The Cingulum: Three ROIs were used in order to recon- 
struct the cingulum bundle. Two seed ROIs were drawn in 
the mid-sagittal plane: first ROI within the body of cingu- 
lum (Figure 1 A, ROI 1), second ROI in the posterior part 
of the cingulum (Figure 1A, ROB). Another was drawn in 
a more lateral sagittal slice in the temporal part of the cin- 
gulum adjacent to the hippocampus (Figure 1A, ROI2). If 
necessary, in order to eliminate fibers which are not part 
of the cingulum, a "no-fiber" plane was drawn in a coronal 
plane posterior to the cingulum and in an axial plane dorsal 
to the cingulum (Figure 1A, no fiber). 

2. The Fornix: two/three ROIs were used to reconstruct the 
fornix. First ROI was drawn on an axial plane, marking 
the body of the fornix (Figure IB, ROI1). Second ROI was 



drawn on the temporal part of the fornix (Figure IB, ROI2). 
If necessary, a third ROI was drawn in the axial plane ante- 
rior to the splenium of the CC. In order to eliminate fibers 
which are not part of the fornix, a "no-fiber" plane was 
drawn in an axial plane dorsal to the fornix (Figure IB, no 
fiber). 

The SLF: two ROIs were used to reconstruct the SLF. First 
ROI was drawn in a sagittal plane in the fronto-parietal part 
of the SLF (Figure 1C, ROI1) and a second ROI was drawn 
in the parieto-temporal part of the SLF (Figure 1C, ROI2). 
The ILF: two ROIs were used to reconstruct the ILF in a sagit- 
tal plane. First ROI was drawn in the occipital (posterior) 
part of the ILF (Figure ID, ROI1). Second ROI was drawn 
in the temporal (anterior) part of the ILF (Figure ID, ROI2). 
If necessary, in order to eliminate fibers which are not part 
of the ILF, a "no-fiber" plane was drawn in a coronal plane 
anterior to the ILF (Figure ID, no fiber). 
The UF: two ROIs were used to reconstruct the UF in a coro- 
nal plane. First ROI was drawn in the temporal part of the 
UF (Figure IE, ROI1), and a second ROI was drawn in the 
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FIGURE 1 | Fiber tracking procedure. Representative seed ROIs drawn on the DTI color coded map, below the plotted fiber bundles for each temporal 
projection. (A) cingulum, (B) fornix, (C) SLR (D) ILF; and (E) uncinate fasciculus. 
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frontal part of the UF (Figure IE, ROI2). If necessary, in 
order to eliminate fibers which are not part of the UF, a "no- 
fiber" plane was drawn in an axial slice dorsal to the frontal 
part of the ILF (Figure IE, no-fiber). 

CREATION OF TRACT ROIs AND REGISTRATION OF DTI MAPS 

1. After the tracking procedure, a mask was created from the 
tracts matrices of all subjects in order to create a tract 
mask. Overall 10 masks were created for each subject: five 
temporal projections (UF, fornix, cingulum, SLF, ILF), in 
both hemispheres, right and left. For each temporal pro- 
jection, the tract masks (that are already normalized to the 
MNI template) of the different subjects were registered to 
a tract mask of one young subject. The registration proce- 
dure was performed using normalized mutual information 
(Studholme et al., 1999) in the SPM software (version 2, 
UCL, London, UK). Spatial registration involved finding 
parameters that either maximize or minimize the normal- 
ized mutual information function, followed by a voxel to 
voxel affine transformation. 

2. The registered tract masks of all subjects were summed to 
create the total tract mask. For each tract only voxels that 
80% of the subjects had a tract passing through them were 
included to create the final fiber ROI. 

3. The same registration parameters that were created for each 
temporal projection were applied (right and left) on the DTI 
maps (ADC, FA). 

4. In order to apply the random field theory (RFT) Spatial 
smoothing with 8 mm full width half maximum Gaussian 
kernel was performed. 

STATISTICAL ANALYSIS 

Voxel-based analysis: VBA is a whole brain technique that allows 
the detection of regionally specific differences in brain tissue com- 
position on a voxel-by-voxel basis (Ashburner and Friston, 2000). 
In this work we used SPM (version 2, UCL, London, UK) to 
perform voxel-by-voxel correlation between the DTI indices and 
performance in each of the cognitive domains. For each tempo- 
ral projection, simple regression was performed between the DTI 
indices (ADC, FA, Da, Dr) and age in each temporal projection 
and multiple regressions were performed between the DTI indices 
(ADC, FA, Da, Dr) co-registered to the appropriate temporal 
projection, and scores in the four cognitive domains (executive 
function, information processing speed, and memory). The mul- 
tiple regressions were controlled for age and motor ability using 
the multiple regression routine in SPM with age and scores in a 
simple finer tapping task were used as constancies. The fiber ROI 
was used as a mask. The statistical threshold was set at p < 0.05, 
after Hochberg sequential correction for multiple comparisons 
(Hochberg, 1988). 

The multiple regression analysis generated statistical para- 
metric maps. The multiple regression routine included applying 
implicit and explicit masks (the appropriate fiber ROI, see above), 
without applying grand mean scaling, threshold masking, or 
global calculation. Correction for multiple comparisons was per- 
formed on the statistical parametric map, and the five ROIs were 
used to constrain the analysis — reducing the multiple comparison 



correction to number of voxels contained within the 10 ROIs 
(5 tracts on L and R sides). The statistical parametric maps are 
presented superimposed on an average Tl image of all subjects, 
with the fiber ROIs drawn upon them allowing an anatomical 
informative reference. 

Scatter plots of scores of the cognitive domains and the aver- 
aged DTI parameters extracted only in the voxels passing the 
significance threshold are presented. A linear fit (y = mx + b) was 
performed and the best linear fit is shown on the scatterplots. 

RESULTS 

Factor analysis of the cognitive tests resulted in three cognitive 
domains: the memory domain, the executive function domain, 
and the information processing speed domain. Significant corre- 
lations between age and cognitive performance were obtained, as 
reported previously (Sasson et al., 2012). 

Positive regression was obtained between age and ADC (p < 
0.05, corrected) in all fibers, including bilateral cingulum, SLF, 
fornix, ILF, and UF. Significant positive regression with age is 
found mostly in frontal, temporal, and parietal parts of the fibers 
and not in occipital parts [Figure 2A Negative regression between 
age and FA (p < 0.05, corrected) was found in all fibers accept 
right SLF, which did not pass the correction for multiple com- 
parisons (Figure 2B)]. Regression between FA and age was more 
prominent in the fornix and in temporal part of the ILF and less 
in the cingulum, SLF and UF. 

Table 2 and Figures 3-7 summarize the results of the mul- 
tiple regression analysis between performance in the cognitive 
domains and DTI parameters corrected for age and motor ability. 
In Figures 3-7, the statistical parametric maps of the partial cor- 
relation between ADC/FA and scores in the cognitive domains in 
each temporal projection corrected for age and motor ability are 
shown (p < 0.05, uncorrected) and below magnification of the 
significant clusters that passed the correction for multiple com- 
parisons (p < 0.05, corrected). In Table 2, only significant voxels 
that passed the statistical threshold (p < 0.05, corrected) are pre- 
sented. Partial correlation between axial and radial diffusivities 
and scores in the cognitive domains are also reported below. 

In Figure 8, scatter plots of scores of the cognitive domains and 
the averaged DTI parameters extracted only in the voxels passing 
the significance threshold are presented. A linear fit (y = mx + b) 
was performed and the best linear fit is shown on the scatterplots. 
In the right cingulum, negative partial correlation was found 
between performance in the memory domain (accuracy scores) 
with ADC in frontal as well as temporal parts of the cingulum 
(Figure 3A) and with Da- In the left cingulum negative partial 
correlation was obtained between performance in the informa- 
tion processing speed domain (RT scores) and FA in parietal part 
of the cingulum (Figure 3B). 

In the right fornix, negative correlation was obtained between 
performance in the memory domain (accuracy scores) with ADC 
(Figure 4A) and with radial diffusivity (Dr), and in bilateral 
fornix between performance in the information processing speed 
(RT scores) domains with FA (Figure 4B). Correlations were 
found mainly in temporal parts of the fornix. 

In the left SLF, negative partial correlation was obtained 
between performance in the executive function domain 
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Table 2 | Partial correlation between DTI parameters: ADC, FA, axial 
diffusivity (Da) and radial diffusivity (Dr) and performance in the 
cognitive domains corrected for age and motor ability (p < 0.05, 
corrected). 



FIGURE 2 | DTI correlates of age. (A) Positive correlation between age 
and ADC (p < 0.05, corrected). Correlations were found in all fibers, 
including bilateral cingulum, SLF, fornix, I LF; and UF Correlation is found in 
frontal, temporal, and parietal parts of the fibers. (B) Negative correlation 
between age and FA (p < 0.05, uncorrected). FA showed less significant 
voxels correlated with age compared to ADC, however, correlation was 
found in all fibers accept right SLF which did not pass the correction for 
multiple comparisons. 



(RT scores) with ADC, Da and Dr and in frontal parts of the SLF 
positive correlation between memory domain (accuracy scores) 
with ADC (Figure 5A) and Dr; in the right SLF negative cor- 
relation was obtained between performance in the information 
processing speed domain (RT scores) and FA in parietal part of 
the SLF (Figure 5B). 

In the left ILF, positive correlation was obtained between 
performance in the memory domain (accuracy scores) and 
FA and positive correlation with Dr in temporal parts of 
the ILF (Figure 6B). In bilateral ILF negative correlation was 
obtained between performance in the information processing 
speed domain (RT scores) and FA in occipital parts of the ILF 
(Figure 6B). 
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Table 2 | Continued 

Executive Information Memory 
function processing 
speed 

ADC 

L 



Uncinate FA 

L 

Dr R 
L 



ILF; inferior longitudinal fasciculus; SLF; superior longitudinal fasciculus; R, right; 
L, left. 

In the right UF negative correlation was obtained between 
performance in the memory domain (accuracy scores) with 
ADC and positive correlation in the executive function 
domain (RT scores) with ADC in temporal parts of the UF 
(Figure 7A). 

In the left UF positive correlation was obtained between 
performance in the information processing speed domain (RT 
scores) with FA in temporal parts of the UF (Figure 7B). 

DISCUSSION AND CONCLUSIONS 

The main finding of the present study is that WM integrity 
is correlated with cognitive performance in a fiber specific 
manner. 

In conclusion the results revealed that performance in the 
executive function domain (RT) is correlated positively with ADC 
and negatively with FA in the left SLF and right UF; performance 
in the information processing speed domain (RT) is correlated 
with negatively with FA in the left cingulum, left fornix, right and 
left ILF and SLF; and the memory domain (accuracy) shows sig- 
nificant positive correlations with FA and negative correlations 
with ADC in the right fornix, right cingulum, left ILF, left SLF, and 
right UF. 

DTI tractography was performed using the free water elimi- 
nation method in order to minimize the partial volume. In the 
correlation analysis we used partial correlation correcting for age 
and motor ability, in order to isolate correlations more specific to 
each of the cognitive domains. 

Functionality has been heavily studied in gray matter regions. 
However, functionality in white matter fiber bundles has only 
recently been investigated. Indeed only in the last 10 years DTI 
tractography enabled us to investigate functionality of white mat- 
ter fiber bundles. Our approach, using factor analysis of a battery 
of cognitive tests, allows us to localize functional domains in 
distinctive white matter fiber tracts. 

Most of our findings are consistent with previous studies and 
a summary of our findings in each temporal projection, as well 
as the known literature on the functionality of the distinctive 



Cingulum 
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FIGURE 3 | DTI correlates of cognitive performance in the cingulum. 

In the upper row, parametric maps of the multiple regression analysis 
between performance in the cognitive domains and DTI parameters (ADC 
and FA) corrected for age and motor ability (p < 0.05, uncorrected) and in 
the lower row, magnification of the significant clusters that passed the 
correction for multiple comparisons (p < 0.05, corrected). (A) Correlation 
with ADC. In the right cingulum, negative correlation was obtained 
between performance in the memory domain with ADC (green clusters). 
(B) Correlation with FA. In the left cingulum correlation was obtained 
between performance in the information processing speed domain and FA 
(yellow clusters). No correlation was obtained between ADC nor FA with 
executive function (red clusters). 



white matter tracts and surrounding regions, as summarized in 
Table 3. 

COGNITIVE PERFORMANCE LOCALIZATION IN THE CINGULUM 

A significant negative correlation was obtained in the right ante- 
rior cingulum between ADC and performance in the memory 
domain. This finding is consistent with previous reports of 
involvement of the anterior cingulate and cingulum in mem- 
ory performance (Valenstein et al., 1987; Sepulcre et al., 2009; 
Charlton et al., 2010). A negative correlation was also obtained 
between ADC in the posterior cingulum with performance in the 
memory domain. 

COGNITIVE PERFORMANCE LOCALIZATION IN THE FORNIX 

Lesion studies suggest that the fornix, which includes choliner- 
gic afferents to the hippocampus, is important for learning and 
formation of new memories (Spiers et al, 2001). 
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FIGURE 4 | DTI correlates of cognitive performance in the fornix. In the 

upper row, parametric maps of the multiple regression analysis between 
performance in the cognitive domains and DTI parameters (ADC and FA) 
corrected for age and motor ability (p < 0.05, uncorrected) and in the lower 
row, magnification of the significant clusters that passed the correction for 
multiple comparisons (p < 0.05, corrected). (A) In the right fornix, negative 
correlation was obtained between performance in the memory domain 
with ADC (green clusters). Positive correlation was obtained in left fornix 
between executive function and ADC, however the clusters did not pass 
the correction for multiple comparisons (red clusters). (B) In bilateral fornix 
between performance in the information processing speed domains with 
FA (yellow clusters). 



Indeed in the present study we obtained that ADC in the right 
fornix is significantly negatively correlated with performance in 
the memory domain. 
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FIGURE 5 | DTI correlates of cognitive performance in the SLF. In the 

upper row, parametric maps of the multiple regression analysis between 
performance in the cognitive domains and DTI parameters (ADC and FA) 
corrected for age and motor ability (p < 0.05, uncorrected) and in the lower 
row, magnification of the significant clusters that passed the correction for 
multiple comparisons (p < 0.05, corrected). (A) Correlation with ADC. In 
the left SLF negative correlation was obtained between performance in the 
executive function domain with ADC (red clusters) and positive correlation 
between memory domain with ADC (green clusters). (B) Correlation with 
FA. In the right SLF negative correlation was obtained between 
performance in the information processing speed domain and FA (yellow 
clusters). 



COGNITIVE PERFORMANCE LOCALIZATION IN THE SUPERIOR 
LONGITUDINAL FASCICULUS 

Various studies suggest an important role for the SLF in higher 
brain functions, and particularly language (Lichtheim, 1885; 
Damasio and Damasio, 1980; Tanabe et al., 1987; Catani et al., 
2003; Wise, 2003; Geldmacher et al., 2007). We obtained a 
negative correlation between performances in the information 
processing speed domain, with bilateral FA in the SLF. The 
correlations were localized to parietal regions within the right 
SLF. The parietal lobe is well known from the literature to be 



related to calculation (Rocha et al., 2005; Cavanna and Trimble, 
2006). 

A strong significant positive correlation was also obtained 
between performance in the executive function domain and ADC 
in the left SLF. The Arcuate fasciculus (AF), a part of the SLF, has 
been central to the neurobiological interpretation of higher brain 
function generally and for language in particular (Lichtheim, 
1885; Geschwind, 1965a,b; Heilman et al, 1970; Tanabe et al, 
1987). The task we used here, stroop interference involves inter- 
ference between semantic information and visual input, and 
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FIGURE 6 | DTI correlates of cognitive performance in the ILF. In the 

upper row, parametric maps of the multiple regression analysis between 
performance in the cognitive domains and DTI parameters (ADC and FA) 
corrected for age and motor ability (p < 0.05, uncorrected) and in the lower 
row, magnification of the significant clusters that passed the correction for 
multiple comparisons (p < 0.05, corrected). (A) Correlation with ADC. In 
the left ILF, positive correlation was obtained between performance in the 
memory domain and FA (green clusters). (B) Correlation with FA. In bilateral 
ILF negative correlation was obtained between performance in the 
information processing speed domain and FA (yellow clusters). No 
correlation was obtained between ADC nor FA with executive function 
(red clusters). 



hence involves linguistic abilities, and is also considered as a 
higher function task. 

COGNITIVE PERFORMANCE LOCALIZATION IN THE INFERIOR 
LONGITUDINAL FASCICULUS 

The ILF is known to play an important role in visual mem- 
ory as demonstrated by post mortem studies (Bauer and Trobe, 
1984) lesion studies (Shinoura et al., 2007) and imaging studies. 
In congenital prosopagnosia correlation was found between face 
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FIGURE 7 | DTI correlates of cognitive performance in the UF. In the 

upper row, parametric maps of the multiple regression analysis between 
performance in the cognitive domains and DTI parameters (ADC and FA) 
corrected for age and motor ability (p < 0.05, uncorrected) and in the lower 
row, magnification of the significant clusters that passed the correction for 
multiple comparisons (p < 0.05, corrected). (A) Correlation with ADC. In 
the right UF negative correlation was obtained between performance in the 
memory domain with ADC (green clusters) and positive correlation in the 
executive function domains with ADC (red clusters). (B) Correlation with 
FA. In the left UF positive correlation was obtained between performance in 
the information processing speed domain and FA (yellow clusters). 



recognition and FA in the ILF (Thomas et al., 2009; Tavor et al., 
2010). In our work we found a significant positive correlation 
between memory performance and FA in the left ILF. The mem- 
ory tasks were visual tasks and included visual spatial memory 
and visual verbal memory. 

There was also negative correlation between FA in the left 
ILF, and performance in the information processing speed 
domain. There are some findings in the literature supporting 
an association between measurements of occipital white matter 
structures with information processing speed. Smaller occipital 
white matter volume was associated with slower information pro- 
cessing speed in schizophrenic patients (Antonova et al., 2005). 
Higher average fractional anisotropy (FA) in left occipitotemporal 
junction was obtained in children trained in mental calculation 
(HuetaL, 2011). 
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FIGURE 8 | Scatter plots of the DTI parameters and scores in the cognitive domains. Scatter plots of DTI parameters extracted from the significant 
clusters (in Y-axis) and Z-scores in the cognitive domains (in the X-axis). 
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Table 3 | A summary of cognitive localization in the temporal projections obtained in the present study along with the known literature. 



Cognitive domains 



Previous findings 



References 



Cingulum Memory 



Fornix 
SLF 



ILF 



UF 



Memory 

Executive function (L, frontal) 
Information processing speed (R, parietal) 



Memory 

Information processing speed 
Memory 



Working memory, attention 
Memory 

Higher brain functions, Language (Uncinate) 
Calculation (parietal lobe) 



Visual memory 



Memory 



Valenstein et al., 1987; Charlton et al., 2006; 
Sepulcre etal., 2009 
Spiers et al., 2001 

Lichtheim, 1885; Geschwind, 1965a,b; 
Heilman et al., 1970; Tanabe et al., 1987; 
Rocha et al., 2005; Cavanna and Trimble, 
2006 

Bauer and Trobe, 1984; Shinoura et al., 
2007; Thomas et al., 2009; Tavor et al., 2010 
Antonova et al., 2005; Hu et al., 2011 
Nakamura etal., 2005; Di Paola etal., 2011; 
O'Dwyer etal., 2011 



COGNITIVE PERFORMANCE LOCALIZATION IN THE UNCINATE 
FASCICULUS 

The UF plays an important role in the formation and retrieval 
of memories as shown by lesion (Levine et al., 1998) and animal 
studies (Squire et al., 2004). Indeed, we found a negative corre- 
lation between ADC in the right UF with performance in the 
memory domain. In a study using diffusion tensor imaging (DTI) 
(using tract based statistics — TBSS) increased axial diffusivity 
was found in the UF in patients with amnestic mild cognitive 
impairment (O'Dwyer et al, 2011). Lower white matter density 
was found in the UF (as well as in the fornix and anterior cin- 
gulum) in patients with hypoxic amnesia (Di Paola et al., 2011). 
Lesion studies have demonstrated that the frontal lobe is cru- 
cial for executive functions (Milner, 1971; Drewe, 1974; Perret, 
1974; Rowe et al, 2001). Indeed, the ADC values in the UF 
which projects into the pre-frontal cortex were obtained in the 
present study to correlate with performance in the executive func- 
tion domain. This finding is consistent with previous imaging 
studies; FA in the UF was correlated with executive function in 
subjects with schizotypal personality disorder (Nakamura et al., 
2005). 

Lower ADC and higher FA imply on higher white matter 
integrity, and indeed in the present paper, they were correlated 
with higher accuracy and lower RT. Higher FA values may indicate 
higher tissue directionality and organization, increased axonal 
and myelin density. Lower ADC values may indicate higher tis- 
sue density (Gouw et al., 2008; Blumenfeld-Katzir et al., 2011). 
Postmortem studies have shown that aging is accompanied by 
loss of 25-45% of the total length of myelinated fibers in the 
brain white matter (Marner et al., 2003). Degeneration of myelin 
sheaths and loss of nerve fibers in several but not all fiber bun- 
dles were correlated with cognitive decline in monkeys, among 
them the fornix (Peters and Kemper, 20 1 1 ). Therefore, it is hard to 
interpret to which biological mechanism to relate the correlation 
of DTI parameters with cognitive decline. Based on the literature, 
it can be attributed to axonal loss, degeneration of myelin sheaths, 
loss of oligodendrocytes (Bartzokis, 2004), or other glial cells. 

Correlation with ADC was accompanied by correlation with 
either Dr, Da or both; however, correlation with FA did not 



correspond with correlation with either diffusivity except in the 
right ILF in which correlation between memory performance and 
Dr was obtained. This can be explained by the fact that ADC is 
an average of both diffusivities; however, change in FA is a result 
of change in the variance between the diffusivities. 

LIMITATIONS OF THE PRESENT STUDY 

DTI used in the present study has several methodological limita- 
tions. Reduced FA and an increase in ADC can result from partial 
volume with the surrounding tissue. In the present study we used 
the free water elimination method (Pasternak et al, 2009) to 
address this issue. The tracking procedure itself has several limita- 
tions. First, crossing fibers cause a reduction in FA and erroneous 
termination of a tract. Second, the fiber tracking procedure is sub- 
jective, and ROI placement can differ among observers. Further, 
imperfect registration of the fiber tracts can be problematic when 
performing the analysis using group voxel-based statistics and 
may include voxels that are not part of the fiber tract. Partial vol- 
ume with surrounding tissue may also result from the smoothing 
procedure. 

SUMMARY 

In conclusion, using DTI tractography, cognitive performance 
can be correlated with WM integrity in the temporal projections. 
Combining voxel wise correlation analysis and fiber tracking, 
enables anatomical definition of ROI for correlation analysis of 
behavioral parameters with diffusion indices, and functionality 
can be correlated with white matter integrity. 
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